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Abstract 



A method for direct sampling of the exponential moments of canonical phase from the 
data recorded in balanced homodyne detection is presented. Analytical expressions for the 
0^ ■ sampling functions are shown which are valid for arbitrary states. A numerical simulation 

illustrates the applicability of the method and compares it with the direct measurement of 
^ ' phase by means of double homodyning. 

§ ■ 1 Introduction 
o 

In the study of the problem of phase of a quantum harmonic oscillator, such as a radiation 
field, one usually proceeds in one of two different ways. In the first, the phase is defined from 



p ' the requirement that phase and photon number should be complementary quantities. This first- 
principle definition leads to the canonical phase related to the one sided unitary phase operator 
[jl], 3. In the second way, phase quantities are defined from the output observed in phase- sensitive 



^ ! measurements, such as eight-port homodyne detection It was found that in such a scheme 
the Q function is measured 0]. The measured phase distribution can then be obtained by radial 
integration of the Q function. Whereas in the classical limit the measured phase coincides with 
^ ' the canonical phase, in the quantum regime the two phases significantly differ from each other in 
■ general. The most important difference is that the integrated Q function yields a broader and less 
structurized distribution than the canonical phase . Is there a way which would bring together 
the advantages of these two approaches - i.e., the theoretical elegance and pronounced structure 
typical of the canonical phase and the advantage of experimental availability as in the case for the 
integrated Q function? 

Before addressing the problem in more detail, let us mention the concept of direct sampling 
of a quantity from experimental data. This concept has been studied extensively in connection 
with the density matrix reconstruction [Q. Assume a balanced homodyne measurement of a field 
quadrature where ^9 is the phase of the local oscillator (LO). Performing such measurements 
with different on large ensembles of identically prepared states we obtain probability distributions 
^9). A quantity A can be sampled from the homodyne data if it can be expressed as a two 
fold integral of the measured distribution, 

A= dd dx, K^{x,^)p{x,^), (1) 

J2n J — 00 



where Ky^^{x,i}) is an integration kernel. The quantity A can represent, e.g., density matrix 
elements, mean values of operators, etc. 

Let us turn to the question of whether the canonical phase distribution can be directly sam- 
pled. In this case the quantity A in Eq. (|l|) is the phase probability, A = p{(p). In Ref. 
it was suggested to obtain the exact phase distribution as the limit of a convergent sequence of 
appropriately parametrized (smeared) distributions each of which can directly be sampled from 
the homodyne data. The exact phase distribution can then be obtained asymptotically to any 
degree of accuracy, if the sequence parameter is chosen such that smearing is suitably weak. How- 
ever, it turned out that whereas the method works well for states with low photon numbers, with 
higher excitation the smearing parameter must be chosen very small and the corresponding ker- 
nels become more and more structurized. This makes sampling problematic for highly excited 
states. We can see a discrepancy between the quantum and classical regions: on one hand phase 
can easily be measured in classical physics, on the other hand sampling of the phase distribution 
becomes tedious for states from the classical region. Can a unified approach be found which would 
bridge the gap between these two regions and which would enable us to measure the canonical 
and classical phase distributions on the same footing? 

In this contribution we show the possibility of direct sampling of the exponential phase mo- 
ments of the canonical phase. The corresponding kernels are well behaved functions which, for 
large \x\, approach their classical counterparts. Since the moments contain the same information 
about the phase properties as the probability distribution itself, the method can serve as a way 
for experimental determination of the canonical phase. 



2 Integration kernels 

The canonical phase distribution p{ip) of a state g is defined as p{ip) = (27r)~^ {(f\g\ip), where 
the phase states \ip) are \ip) = Z].^o ^^"'^1'^)' l'^) being the Fock states. The exponential phase 
moments ^'/j of this distribution are given by = e"^"^ p{ip) dip and they can be expressed as 

oo 
n=0 

for k positive and = "^-k fo^^ ^ negative. Our aim is to express by means of the measured 
quadrature distribution p{x,'d), i.e., in the form of Eq. (|I|) with A = "^k- For this purpose we 
must find the corresponding kernel Kk[x,i)). 

To do so, let us express the distribution p{x, {}) by means of the density matrix elements as 

oo oo 
rp=0 nv=0 

Here i/^nix) are the eigenfunctions of the harmonic oscillator Hamiltonian, ipni^) = {2""n\^y^l)^^^^ 
exp(— a;^/2) H„(x), H„(x) being the Hermite polynomials. Substituting Eq. (D into (|1|) and 
comparing with Eq. (|^) we find that the kernels must be of the form Kk{x, = e*^'' Kk{x), where 
the x-dependent part must satisfy the integral equation 

/oo 
dx Kk{x)^p^n^{x)ipnix) = 1 (4) 
-oo 



(n = 0, 1, 2, . . .). In P the solution of this equation is discussed in detail Let us write here the 
result in the form 

K2m{x) = -, ^ TT / ati e 1 ■ ■ ■ 
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with Zfc = [exp(-Ej=itj^) 
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l]/2, ^{a,b,y) being the confluent hypergeometric function. 
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Figure 1: The functions Kk{x). 

In Fig. 1 we plot the functions Kk{x) for several values of k. As can be seen, the functions 
Kk{x) are well behaved and with increasing \x\ they quickly approach their (classical) asymptotics 
K^{x), where 



and 
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being a (unimportant) constant. An essential difference between Kk{x) and K^{x) appears 
only for x near zero, within the area of vacuum fluctuations. It is shown in p] that the functions 
K^{x) can be obtained as kernels for sampling exponential phase moments in classical physics. 
Therefore, we have found kernels for direct sampling of the canonical phase moments which can 
be used for any state, regardless if it has typically quantum or classical properties. 



3 Simulated measurements 



To illustrate the applicability of the method we have performed computer simulations of homodyne 
measurements and used the kernels Kk{x,d) to determine the moments '^k- For comparison, 
we have also simulated double homodyne measurements to get exponential phase moments that 
correspond to the radially integrated Q function. In the computer simulations the state to be 
detected is the phase squeezed state |a, s) with the coherent amplitude a = 5xe^'^°, ipo = 0.6, and 
the squeeze parameter s = 6 (the mean photon number of this state is (n) = 26.04). For both the 
homodyne and the double homodyne measurements the total number of measurement events is A^^e 
= 6020. For the homodyne measurements the LO phase takes 41 values equidistantly distributed 
over the 27r interval. As shown in , the statistical error of the sampled phase moments depends 
on the numbers of measurement events for different i). By a proper distribution of the total number 
iVe for individual phases we can decrease the statistical error of various moments In order 
to minimize the statistical error of the first moment we increased the number of measurement 
events for such -(9 for which the peak of p{x^d) is near x = (maximum 800 events), whereas for 
d yielding a peak far from zero the number of events was small (minimum 10 events). 
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Figure 2: Real (a) and imaginary (b) parts of the experimentally determined exponential phase 
moments Bars with full lines: direct sampling from balanced homodyning, bars with dashed 
lines: double homodyning. The vertical lines represent the estimated statistical error. 



The experimentally determined phase moments together with the estimated statistical errors 
are shown in Fig. 2. We can see that the absolute values of moments of the integrated Q function 
are smaller than the corresponding values of the canonical distribution; the difference becomes 
larger with increasing k. This corresponds to the fact that the integrated Q function smears the 
structure of the canonical distribution. 

Let us compare the statistical errors. For the odd moments the errors are approximately of 
the same magnitude in the two methods, whereas for the even moments the sampling method 
yields larger errors. This reflects the qualitatively different behavior of the kernels for k odd and 
k even. It is also related to the chosen numbers of measurement events for different d: in our 



example we have distributed the event numbers so as to minimize which on the other hand 
increases the errors of even moments. Even though one could expect that the double homodyning 
- as a direct phase measurement - would yield generally smaller statistical errors, we find that the 
statistical error of the first moment is smaller for the sampled canonical distribution. (The 
direct measurement means that a single measurement event yields a single value of phase.) 

Let us mention that the first moment is connected to very important characteristics of the 
phase distribution. The mean value of phase (p can be calculated d& (p = arg\l/i; this quantity can 
correspond, e.g., to a phase shift in an interferometer. Since is determined more precisely in the 
balanced homodyning than in the double homodyning, the sampling method enables us to deter- 
mine ^ with smaller statistical error. From the experimental data we obtain = 0.5994±0.0011 
for the sampling method, whereas from the integrated Q function we obtain = 0.5990±0.0016. 
(Note that for both distributions the correct value is = y^o = 0.6.) As can be seen, the error 
of determination of the mean phase by means of homodyne sampling is about 70% of the error 
in the double homodyning. The moment is also related to various phase uncertainties, which 
describe the "width" of the phase probability distribution. A phase uncertainty Aty? can be defined 
as A(/9 = arccos|\I^i|, which is related to the Bandilla-Paul phase dispersion a\p as a bp = sinAy9 
and to the Holevo phase dispersion af^ as an = tanAip (An advantage of the uncertainty Aip 
is that it enables us to measure the phase width in the same units as the phase itself - in radians, 
degrees, etc.) In this way we obtain Aip = 0.065 for the canonical phase distribution and Aip = 
0.125 for the integrated Q function. 

4 Discussion and conclusion 

The presented method shows a very straightforward way for obtaining the exponential moments 
of the canonical phase from the data of homodyne detection. Direct sampling enables us to 
reconstruct the moments in real time as the experiment runs, together with the estimation of 
the statistical error |^ . In this way the theoretically profound concept of canonical phase can be 
connected with data obtained from present experiments. 

The moments contain the same information as the probability distribution p{<p). Therefore 
they can be used for reconstruction of the original function pi^p). However, even the lowest 
moments give us an interesting information about the phase properties, e.g., the first moment 
is directly related to the mean value of phase and to the phase uncertainty. 

The integration kernels are well-behaved functions which rapidly approach their asymptotics 
given either as step- functions (for odd moments) or logarithmic functions (for even moments). As 
shown in [Q, these functions can serve as kernels for sampling of the phase moments in classical 
physics. Therefore, we have found a unified approach which connects the measurement of the 
canonical phase with its classical counterpart. 

It has been shown that the accuracy of the sampled moments of the canonical phase is com- 
parable with that of the directly measured radially integrated Q function. Moreover, when the 
total number of measurement events is the same, the sampling method can yield the mean value 
of phase more precisely than the measurement of the Q function. Also this aspect can make the 
presented method very attractive for experimental applications. 
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